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PREFACE 


In  this  paper,  we  describe  the  results  of  some  computerized 
Monte  Carlo  simulations  of  a stochastic  attrition  process  anal- 
ogous to  the  Lanchester  "linear-law"  differential-equation  model 
of  combat.  The  problems  studied  include  expected  numbers  and 
distributions  of  survivors  as  functions  of  time,  the  effects  of 
iterative  approximations  to  expected  numbers  of  survivors,  and 
the  probability  of  each  side's  "winning"  an  engagement  in  which 
the  termination  rule  is  of  the  form  of  breakpoints  (not  neces- 
sarily equal)  for  each  side.  In  the  latter  case,  we  also  con- 
sider distributions  of  numbers  of  survivors. 
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Chapter  I 

THE  MATHEMATICAL  MODEL 


The  deterministic  model  of  which  the  stochastic  process 
studied  here  is  an  analogue  was  introduced  by  F.  W.  Lanchester 
(Reference  [6]).  After  presenting  his  "square  law"  as  a model 
of  combat  in  which  the  numerically  superior  side  is  able  to 
bring  that  superiority  to  bear  on  the  opposition,  Lanchester 
turned  to  the  description  of  combat  that  occurs  in  the  sense 
of  one-on-one  engagements,  so  that  the  numerically  superior 
side  has  an  advantage  only  in  having  more  eligible  combatant  s._j 
Lanchester  proposed  the  model 

b"(t)  = - c^CtJrCt)  and  r"(t)  = - c2b(t)r(t)  , (1) 

where  b(t)  and  r(t)  are  the  numbers  of  surviving  Blue  and  Red 
combatants,  respectively,  t time  units  after  the  combat  begins. 
Here  c^  and  c2  are  positive  constants. 

We  make  the  following  observations  for  purposes  of  analysis 
of  simulation  results: 

If 

c-jTCO)  = c2b(0)  , 

then 

lim  b ( t ) = lim  r(t)  = 0. 
t-*00  t-*00 

Otherwise,  lim  b(t)  > 0 and  lim  r(t)  = 0 or  lim  b(t)  = 0 and 
lim  r (t)  > 0 according  as  c^r(0)  < c2b(0)  or  c^r(0)  > c2b(0). 
Thus,  the  force  ratio 
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0^(0) 
f = c^bTO) 

determines  the  victor  when  the  engagement  is  continued  until  one 
side  or  the  other  is  annihilated. 

If 


C1  c2  * 

then 

b(t)  - r (t ) = b(0)  - r (0)  (2) 

for  all  t.  Equation  (2)  confirms  the  nature  of  the  model  as  a 
description  of  combat  occurring  in  the  form  of  one-on-one  duels 
and  of  the  role  of  numerical  superiority.  Here  the  originally 
superior  side  is  unable  to  increase  its  absolute  superiority, 
which  it  is  able  to  do  in  square-law  combat;  cf.  Reference  [5, 

p.  10]. 

The  stochastic  process  we  have  simulated  is  the  homogeneous 
linear-law  process  LI  of  Reference  [4].  For  completeness,  we  list 
here  the  assumptions  from  which  that  process  is  derived,  along 
with  a description  of  the  process  as  a regular  Markov  process. 


Assumptions 

(1)  All  combatants  on  each  side  are  identical. 

(2)  The  time  required  for  a particular  Blue  combatant  to 
detect  a particular  Red  combatant  is  exponentially 
distributed  with  expectation  1/d^j  Each  Blue  combat- 
ant detects  different  Red  combatants  independently. 

(3)  A Blue  combatant  attacks  every  Red  combatant  it  detects; 
the  conditional  probability  of  kill,  given  detection  and 
attack,  is  k^ . The  attack  occurs  instantaneously,  and 

contact  is  lost  immediately  thereafter.  No  attack  can 
occur  without  a detection. 

(4)  Red  combatants  satisfy  the  same  assumptions  with  param- 
eters d„  and  k . 

r r 

(5)  The  detection  and  attack  processes  of  all  combatants 
are  mutually  independent. 


2 


One  can  construct  a probability  space  (ft,M,P)  and  on  it  a 
stochastic  process  ( ( B^. , ^ t > q state  space  E = N x N 

(where  N = {0,1, 2, •••,})  such  that  the  following  characteriza- 
tion holds.  The  interpretation  is  that  Bt  is  the  (now  random) 
number  of  Blue  combatants  surviving  at  time  t and  that  R.  is  the 

U 

.corresponding  number  of  Red  combatants. 

THEOREM.  Under  Assumptions  (l)-(5)>  the  stochastic  process 
( (B^ ,Rfc ) )t >Q  is  a regular  Markov  process,  with  jump  function  X 
given  by 

X(i,j)  = ij(Vb+drkr)  , 
transition  kernel  P given  by 

P((i,j>i(i-i.J))  - d ,Xrd  k 

b b r r 

P((i>J  ) i (i>j-l)}  = h k + d k * ^3) 

b b r r 

and  infinitesimal  generator  Q given  by 

Q (( i » J ) J ( i~l , J ) ) = iJdrkr 

Q((i,j);(i,j))  = -ij(dbkb+drkr) 

Q((i»j ) ; (i,j-l))  = iJdbkb  . (4) 

The  assertion  that  this  stochastic  process  is  an  appropriate 
analogue  of  the  deterministic  model  (1)  is  based  upon  the 
resemblance  of  (1)  and  the  infinitesimal  generator  Q given  in 
Equation  (4).  (Reference  [4]  gives  a detailed  explanation  of 
this  analogy. ) 

The  reader  is  referred  to  Reference  [3]  for  further  details, 
interpretations,  and  proofs  concerning  regular  Markov  processes. 
By  a well-known  characterization  of  regular  Markov  processes 
(Reference  [1]),  we  have  the  following  interpretation:  When 

the  stochastic  attrition  process  ((Bt,Rt))t>Q  enters  a state 
(i,j),  it  remains  there  an  exponentially  distributed  sojourn 
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which  is  independent  of  all  past  history  of  the  process  and  has 
expectation  1/A(i,j).  At  the  end  of  this  sojourn,  it  jumps  to 
a new  state  according  to  the  distribution  P((i,J);.)  indepen- 
dently of  all  past  history.  Including  the  length  of  the  sojourn. 

A A 

The  sequence  ((Bn,Rn))n^  of  states  visited  forms  a Markov  chain 
(the  imbedded  Markov  chain)  with  transition  kernel  P given  by 
Equations  (3).  Except  along  the  axes,  this  Markov  chain  is  a 
homogeneous  random  walk — an  observation  that  is  useful  in  the 
analysis  of  some  of  our  simulation  results. 

We  denote  by  P^’^  the  probability  law  of  the  attrition 
process  ((B^,!^))  subject  to  the  initial  conditions 

BQ  = i and  RQ  = j 

and  by  the  corresponding  expectation  operator. 

The  details  of  the  simulation  procedure  are  essentially 
those  used  in  a previous  simulation  study  of  a stochastic 
square-law  attrition  process  (we  refer  the  reader  to  Reference 
[5]  for  details). 
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Chapter  II 

EXPECTED  NUMBERS  AND  DISTRIBUTIONS  OF  SURVIVORS 


We  begin  with  a rather  detailed  attempt  to  ascertain  the 
shapes  of  the  functions 

t -v  E(i,J)[Bt]  and  t -*•  E(1,J)[Rt] 

for  various  choices  of  initial  force  strengths  (i,j)  and  fixed 

values  of  the  parameters  d,  , d , , and  k . In  Table  1,  we 

b r*  b r 

summarize  the  results  obtained. 

Table  1 contains  three  sets  of  data,  corresponding  to 


Bo 

= 50 

and 

II 

o 

cc 

60  ; 

Bo 

= 50 

and 

Ro  = 

55  . 

The  most  striking  aspect  of  all  these  results  is  the  great 
intensity  of  the  early  stages  of  the  engagement.  In  Figure  1, 
we  give  graphically  the  progress  of  the  battle.  Since,  for  all 
t,  we  very  nearly  have  E[Bfe]  = E[Rt]  (the  two  forces  are  ini- 
tially exactly  equal),  we  have  given  one  graph  for  both  sets  of 
data.  One  can  observe  that  each  side  is  reduced  to  half  its 
initial  strength  by  t = 0.1,  even  though  the  mean  time  to 
detect  a particular  member  of  the  opposition  is  1/d^  = 1/d^  * 
10.  Initial  intensity  is  greater  for  linear-law  combat  than 
for  square-law  combat  with  appropriately  analogous  parameters. 

Figure  2 represents  the  same  data,  but  with  the  time  axis 
logarithmic.  The  figure  shows  that  the  attrition  begins  very 
intensely,  and  then  moderates  as  the  engagement  progresses. 
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Table  1.  EXPECTED  NUMBERS  OF  SURVIVORS 


[d 


b 


d =0.1 
r 


k =0.5] 
r J 


t 

E(l 

E ^ 1 ,j)[Rt] 

t 

E(i'J)[Bt] 

E(l  ,J )[Rt] 

0.000 

0.100 

0.200 

0.300 

0.400 

0.500 

0.600 

0.000 

0.010 

0.020 

0.030 

0.040 

0.050 

0.060 

0.070 

0.080 

0.090 

0.100 

0.000 

0.001 

0.002 

0.003 

0.004 

0.005 

0.006 

0.007 

0.008 

0.009 

0.010 

0.000 
1 .000 
2.000 

3.000 

4.000 

0.000 

0.100 

0.200 

0.300 

0.400 

0.500 

0.000 

0.010 

0.020 

0.030 

200.00 

99.20 

67.36 

49.48 
4 0.34 
34.22 
29.86 

200.00 
181  .40 

166.96 
153.44 
143.12 
1 33.02 
124.66 
117.74 
111.18 
105.34. 
100.10 

200.00 

197.96  , 
195.76 

194.58 
192.48 
1 90.50 
188.40 
187.42 

184.96 
183.04 

182.58 

50.00 

10.60 

4.96 

3.18 

2.44 

50. DO- 
38. 56 
30.02 
25.18 
22.24 

18.48 

50.00 
48.32 
4 7.08 
45.98 

200.00 

100.64 
65.94 
52.64 
41  .20 
33.52 
30.30 

200.00 
181 .86 
165.98 
153.22 
142.66 
133.62 
125.80 
118.00 

110.64 
106.44 
100.10 

200.00 

198.28 

196.28 
194.56 
192.34 
190.26 

189.02 

187.02 
186.10 
183.86 
181.52 

60.00 
19.86 
14.72 
1 3.04 
11.62 

60.00 

48.18 
41  .52 
35.84 
30.14 
28.10 

60.00 

58.58 

57.18 
56.04 

0.040 

0.050 

0.060 

0.070 

0.080 

0.090 

0.100 

0.000 
0.001 
0.002 
0 . 003 
0.004 
0.005 
0.006 
0.007 
0.008 
0.009 
0.010 

0.000 
1 .000 
2.000 

3.000 

4.000 

5.000 

0.000 

0.100 

0.200 

0.300 

0.400 

0.500 

0.600 

0.700 

0.800 

0.900 

0.000 

0.010 

0.020 

0.030 

0.040 

0.050 

0.060 

0.070 

0.080 

0.090 

0.100 

44.78 

43.58 
42.54 
41  .20 
39.98 

39.20 

38.76 

50.00 

49.84 

49.76 
49.64 

49.40 
49.34 
49.18 
48.82 
48.80 

48.84 
48.72 

50.00 

13.04 

6.28 

4.70 

3.54 

3.16 

50.00 
38.86 
31  .76 
26.52 
23.02 
20.68 
17.36 

16.58 

14.94 

13.78 

50.00 

48.72 

47.28 

46.96 

44.90 

44.24 

42.76 
41  .86 

40.40 

40.20 

39.94 

55.22 

53.02 
51  .72 
51  .74 
50.20 
49.42 

47.54 

60.00 

59.80 

59.68 

59.62 

59.44 

59.36 

59.16 
58.88 
58.74 
58.76 

58.58 

55.00 

16.86 

12.92 

10.14 
8.36 
8.08 

55.00 

44.26 
36.34 

32.14 
28.24 
25.32 

23.92 
22.04 

19.58 

18.22 

55.00 

54.02 

52.26 
51  .42 

50.14 

49.16 
48.10 
46.66 

46.54 
44.86 

43.26 
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Figure  1.  EXPECTED  NUMBERS  OF  SURVIVORS  (DETECTION  RATE  = 0.1) 
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Figure  2.  EXPECTED  NUMBERS  OF  SURVIVORS  (Logarithmic  Scale) 


Neither  graph  Indicates  a simple  form  for  the  function  t -*• 

E[B,]  = E[R. ].  It  is  not,  for  example,  a linear  function  of 

1/  U 

log  t.  One  reason  for  undertaking  these  investigations  that  use 
Monte  Carlo  simulations  was  our  inability  to  obtain  a tractable 
analytical  expression  for  these  functions  (cf.  Reference  [5]  for 
some  pertinent  comments  concerning  the  difficulties  involved). 
The  simulations  have  given  us  a qualitatively  clear  picture 
(intense  attrition  at  the  start,  becoming  less  intense  as  time 
elapses),  but  the  analytical  problems  remain  and  are  worthy  of 
further  research  effort. 

The  other  two  sets  of  data  in  Table  1 exhibit  substantially 
the  same  qualitative  behavior  and  in  this  respect  will  not  be 
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discussed  further.  A more  interesting  aspect  of  these  data  is 
that  they  seem  to  exhibit  the  stochastic  analogue  (in  the  sense 
of  expectations)  of  the  property,  expressed  in  Equation  (2), 
of  preservation  of  absolute  numerical  advantage.  Namely,  in 
both  cases,  we  nearly  have 

E(1,J)rBt  - V - Bo  - Ro  - 1 - J 

for  all  t.  (See  Figure  3 for  a graphical  presentation  of  the 
third  set  of  data  in  Table  1.)  We  are  led,  therefore,  to  the 
following: 

CONJECTURE.  If  d.k.  * d k , then  for  all  initial  conditions 

b b r r 

(i,j)  and  all  t > 0 we  have 

E(i,J)[Bt  - Rt]  - i - j . (5) 

To  see  why  this  conjecture  should  be  true,  consider  the 

/S  A 

imbedded  Markov  chain  ((Bn,Rn)).  We  warn  the  reader  that 

(B  »R  ) / ( B » R ) > indeed,  (6n>Rn)  * (BT  ,RT  ),  where  Tn  is 

t h n n 

the  random  time  of  the  nu  change  of  state  of  the  continuous 
time  process  ((B  ,Rt  ))  . In  this  case,  by  Equations  (3), 

P((i,j);(i-l,j))  = P((i,j  ) = (i, j-D)  = 1/2 

shows  that  we  have,  with  respect  to  P^*^,  the  representation 

(Bn,Rn)  = (i,j)  + | Uk,Yk),  n < min{iTj}  , 

where  (X1,Y1)  and  (X,,,Y2),*’*  are  independent  and  identically 
distributed  random  variables  with 

p(i,j  ){(Xk>Yk)  = (-i,0)}  = P(1  J){(Xk,Yk)  = (0,-l)}  = 1/2  . 
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Figure  3.  EXPECTED  NUMBERS  OF  SURVIVORS  (Unequal  Forces) 
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Thus,  it  is  clear  that 

E(1-J)[Bn  - Rn]  = i - j (6) 

for  n <_  min{i,j).  For  other  n,  the  argument  becomes  more  com- 
plicated, but  Equation  (6)  remains  true.  The  only  step  remain- 
ing is  to  pass  from  Equation  (6)  to  Equation  (5),  which  should 
not  be  too  difficult. 

In  Figures  4,  5,  and  6,  we  give  the  probability  distribu- 
tion of  the  state  of  the  process  at  times  0.3125,  1.25,  and  5 
for  the  parameter  values 


a particular  location 
the  simulation  for  which 


results  seem  warranted: 

= 0.3125,  in  Figure  4), 
in  the  vector  random 
(Bt,R^).  There  is  also  rather  evident  symmetry  in 

its  distribution.  As  time  increases,  the  variance,  as 
one  would  certainly  predict,  increases  also. 

(2)  For  t = 0.3125  (Figure  4)  and  t = 1.25  (Figure  5), 
there  is  a negative  correlation  between  Bt  and  R , 

which  is  (visually,  at  least)  remarkably  close  to 
-1.  One  might  have  so  expected  on  intuitive  grounds; 
nonetheless,  it  is  illuminating  to  have  a confirmation. 
For  t = 5 (Figure  6),  the  negative  correlation  is 
present  but  not  -1,  because  of  the  possibility  (indeed, 
frequency)  of  annihilation.  To  us,  this  negative  cor- 
relation is  a rather  strong  argument  against  the  use  of 
deterministic  models,  even  though  (as  seen  in  the  next 
chapter)  the  symmetry  of  these  distributions  allows 
iterative  calculations  of  expected  numbers  of  survivors. 
This  is  because  (with  initially  equal  forces)  in  a 
deterministic  model  attritions  to  the  two  sides  are 
always  equal;  that  is,  have  correlation  equal  to  1. 


B0  = Rq  = 100 

db  - dr  = 0-1 

k.  = k = 0.5 
b r 

In  all  figures,  the  digit  appearing  at 
(x,y)  is  the  number  of  realizations  of 
(Bfc,Rt)  = (x,y). 

Two  observations  concerning  these 

(1)  Even  for  small  times  (i.e.,  t 
there  is  substantial  variance 
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Clearly  negatively  correlated  attritions  are  more 
physically  plausible. 

The  reader  may  draw  his  own  further  conclusions  from 
these  data. 
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Figure  4.  PROBABILITY  DISTRIBUTION  OF  NUMBERS  OF  SURVIVORS 
(Time  = 0.3125) 
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Chapter  III 
ITERATED  CALCULATIONS 

As  discussed  in  more  detail  in  Reference  [5],  it  is  of 
some  interest  in  the  context  of  the  potential  applicability 
of  such  stochastic  attrition  models  (through  the  use  of  ex- 
pectations) to  iterative  and  deterministic  combat  simulations 
(e.g.,  IDAGAM  I — Reference  [1])  to  consider  the  error  in  the 
approximation 

E(i,J)[Btl+t2]  = E(b ^ [B^  ] , (7) 

where 

(b,r)  = E(i’J)[(B.  ,R.  )] 

tx  z1 

= (E(i5j)[B.  ],  E(i,J)[R.  ])  . 

bl  bl 

In  order  that  the  question  even  make  sense,  one  must  extend 
the  state  space  of  the  stochastic  attrition  process  to  be 
[0,°°)  x [0,«°).  This  extension  represents  no  real  problem 
(details  are  the  same  as  described  in  Reference  [5]  for  the 
square-law  process). 

The  results  of  these  investigations,  for  several  choices 
of  parameter  values  and  initial  forces,  appear  in  Table  2. 

The  five  sets  of  data  exhibit  the  same  qualitative  features, 
the  three  most  salient  of  which  are  the  following: 

(1)  In  absolute  terms,  the  error  committed  by  using  the 
iterated  calculation  to  approximate  the  correct  cal- 
culation is  small — never  more  than  about  3*  In 
relative  terms,  the  error  may  be  substantial 
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Table  2.  ITERATED  CALCULATIONS  OF  NUMBERS  OF  SURVIVORS 


B 

= R 

o Ko 

100;  db 

= dr  * i 

0.1;  kb 

* kp  « 0 

.5 

Iteration 

Interval 

0. 

5 

1 . 

0 

2. 

0 

4. 

0 

Time 

Blue 

Red 

Blue 

Red 

Blue 

Red 

Blue 

Red 

0.0 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

100.00 

0.5 

28.08 

28.88 

-- 

-- 

-- 

-- 

-- 

-- 

1.0 

16.62 

16.82 

17.54 

17.32 

-- 

-- 

-- 

1 .5 

11  .70 

11  .62 

-- 

-- 

-- 

-- 

-- 

-- 

2.0 

9.42 

8.80 

9.86 

9.02 

10.66 

11.10 

-- 

-- 

2.5 

7.78 

7.10 

-- 

-- 

-- 

-- 

-- 

-- 

3.0 

6.62 

5.96 

7.02 

5.82 

-- 

-- 

-- 

-- 

3.5 

5.90 

5.22 

-- 

-- 

-- 

-- 

-- 

-- 

4.0 

5.10 

4.84 

5.36 

4.58 

5.14 

5.52 

7.82 

6.48 

Bo  - 

100;  Rq 

« 125; 

d.  = d 
b r 

* 0.1;  k 

* k * 
b r 

0.5 

0 

- . 

- . 

100.00 

125.00 

100.00 

125.00 

100.00 

125.00 

1 

8.06 

33.96 

-- 

-- 

-- 

-- 

2 

-- 

-- 

2.14 

27.20 

31  .60 

26.10 

-- 

-- 

3 

0.84 

25.56 

-- 

4 

-- 

-- 

0.18 

25.06 

0.26 

23.18 

0.68 

26.26 

Bq  - 100;  Rq  = 

125;  db 

= d = 
r 

0.1;  kb 

* 0.5;  k 

r * 0.4 

0 

100.00 

1 25.00 

100.00 

125.00 

100.00 

125.00 

1 

-- 

16.42 

21.94 

-- 

-- 

-- 

-- 

2 

8.60 

12.18 

10.50 

13.10 

-- 

-- 

3 

-- 

5.50 

9.18 

-- 

-- 

-- 

-- 

4 

-- 

-- 

4.10 

6.76 

5.18 

6.28 

6.06 

10.70 

B0  ' 

100;  Rq 

« 150;  d 

b = dr  " 

0.1;  kb 

* 0.5; 

kp  * 0.33 

0 

-- 

-- 

100.00 

150.00 

100.00 

150.00 

100.00 

150.00 

1 

-- 

-- 

16.48 

27.42 

-- 

-- 

-- 

-- 

2 

-- 

-- 

8.24 

16.08 

9.62 

15.92 

-- 

-- 

3 

-- 

-- 

5.34 

11  .50 

-- 

-- 

-- 

-- 

4 

-- 

-- 

3.82 

9.10 

4.42 

8.50 

6.50 

10.82 

Bo  - 

200;  Rq 

, - 250; 

db  = dr 

= 0.1;  k 

b - °-5- 

kr  = 0 

.4 

0 

-- 

-- 

200.00 

250.00 

200.00 

250.00 

200.00 

250.00 

1 

-- 

-- 

17.67 

27.78 

-- 

-- 

-- 

-- 

2 

-- 

8.20 

15.42 

12.66 

14.04 

-- 

-- 

3 

-- 

4.88 

11  .26 

-- 

-- 

-- 

-- 

4 

-- 

-- 

3.16 

9.24 

6.30 

6.12 

10.58 

11  .06 
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compared  to  surviving  forces,  but  is  never  large 
relative  to  initial  numbers  of  forces. 

(2)  The  iterated  calculation  always  underestimates  the 
value  obtained  in  the  proper  calculation,  and  the 
amount  of  underestimation  increases  with  the  number 
of  iterations. 

(3)  The  amount  of  underestimation  by  iterative  calcula- 
tions increases  with  time.  Hence,  for  times  at 
which  attrition  is  small  (compared  to  initial  forces), 
the  approximation  is  probably  accurate  enough  for 
most  applications.  A difficulty,  though,  is  that 
only  over  very  short  time  periods  is  the  relative 
attrition  fairly  low,  because  of  the  great  initial 
intensity  of  linear  law  combat. 

We  have  no  theoretical  explanation  of  the  underestimates 
arising  through  iterative  calculations. 
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Chapter  IV 

BREAKPOINTS  AND  FORCE  EQUALITY 


Here  we  consider  the  following  situation:  Suppose  that 

for  Blue  and  Red  we  have  defined  breakpoints  and  q^  (each 

between  zero  and  one)  and  that  the  combat  terminates  whenever 

IV  < q,B_  or  R,  < q R_.  Further,  we  declare  Red  to  be  the 
t — b 0 t — r 0 

winner  if  Bfc  <_  qbBQ  occurs  first.  For  various  choices  of 
force  parameters  and  breakpoints,  we  are  interested,  then,  in 
computing  the  probability  of  Red’s  winning. 

In  analytical  terms,  we  seek  to  compute 

p(i, j ){U  < v)  , 

where  U and  V are  the  stopping  times  defined  by 

U = inf  (t : Bt  < qbBQ}  and  V = inf(t:  R < qrRQ}  . 

This  probability,  however,  can  clearly  be  computed  by  using  the 
imbedded  Markov  chain  ( ( Bn , Rn ) ) n> q 3 terms  which  it  is 

expressed  as 

p(i, j ) {s  < T}  , 

where 

A A A A 

S = inf{n:  B < q.B  } and.  T = inf{n:  R'  < q R_)  . 

n — b 0 n — r 0 

Because  of  the  particularly  simple  structure  of  the  imbedded 
Markov  chain,  one  should  be  able  to  compute  these  probabilities 
explicitly — a problem  on  which  we  plan  future  research. 

We  have  obtained  simulation  results  for  cases  of  equal  and 
unequal  breakpoints,  in  which  we  fix  force  numbers  and  compute 
the  probability  of  Red’s  winning  as  a function  of  the  Lanohester- 
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linear  force  ratio 


f 


drkrR0 

Wo  ’ 


(8) 


which  is  varied  by  changing  the  value  of  the  Red  detection  time 
1/d^.  The  Lanchester-linear  force  ratio  is,  of  course,  the 
same  quantity  usually  called  just  the  "force  ratio."  We  dis- 
tinguish it  because  the  appropriate  force  ratio  depends  on  the 
stochastic  model  under  consideration  (for  details,  cf. 

Reference  [5] ) • 


The  first  breakpoints  considered  are  = q^  = 0.5  (the 
results  are  given  in  Table  3 below;  some  of  these  also  appear 
in  Figure  7 ) . 


The  main  conclusions  to  be  drawn  from  these  results  are 
the  following: 


(1)  For  i = j , as  i increases  the  graph  of  the  function 

f ) (Red  wins}  = p^(f)  becomes  different  from  0 

or  1 over  a decreasingly  small  set  that  contains  the 
point  1 (at  which  we  clearly  have  each  side  equally 
likely  to  win).  Indeed,  it  is  easy  to  show  that 


lim  p (f ) 

i-voo  1 


0, 

1/2, 

1, 


f < 1/2 
f = 1/2 
f > 1/2  . 


The  graph  is  also,  consequently,  increasingly  steep 
over  the  set  where  it  is  neither  zero  nor  one. 


(2)  The  force  ratio  required  to  achieve  a substantial  (say 
0.9)  probability  of  victory  is  much  less  than  that 
postulated  by  the  "conventional  wisdom,"  which  holds 
that  a force  ratio  of  2 or  3 is  required  for  victory. 
In  all  the  cases  shown  in  Figure  7,  a force  ratio  of 
1.4  gives  a 0.9  probability  of  victory.  The  larger 
the  number  of  combatants  (on  both  sides)  the  smaller 
the  force  ratio  required  for  0.9  probability  of 
victory.  These  conclusions  are  of  significant 
practical  import. 


In  Table  4 and  Figure  8,  we  present  analogous  results  for 
the  asymmetric  breakpoints 
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Table  3.  PROBABILITY  OF  WIN  AS  A FUNCTION  OF 
FORCE  RATIO  (Equal  Breakpoints) 

[qb  = qr  = 0-5] 


f 

P ( Red  Wins) 

f 

P(Red  Wins) 

f 

P ( Red  Wins) 

Bl 

■ Ro  ■ 50 

Bo 

= R0  = 75 

Bo  ' 

: 50 

’ R0  = 75 

0.5 

0.01 

0.5 

0.01 

0.240 

0-00 

0.6 

0.02 

0.6 

0.01 

0.489 

0.00 

0.7 

0.11 

0.7 

0.07 

0.738 

0.15 

0.8 

0.18 

0.8 

0.19 

0.987 

0.48 

0.9 

0.32 

0.9 

0.35 

1 .236 

0.81 

1 .0 

0.47 

1 .0 

0.50 

1 .485 

0.97 

1 .1 

0.65 

1.1 

0.71 

1.734 

0.99 

1 .2 

0.72 

1 .2 

0.80 

1 .983 

1 .00 

1 .3 

0.85 

1.3 

0.86 

-- 

1 .4 

0.84 

1.4 

0.92 

- - 

Bo' 

50,  Rq  = 100 

Bo  ■ 

Rq  = 200 

Bo  = 

Ro 

= 500 

0.4 

0 -00 

0.75 

0.01 

0.75 

0 .00 

0.5 

0.00 

0.80 

0.05 

0.80 

0.00 

0.6 

0.03 

0.85 

0.16 

0.85 

0.04 

0.7 

0.09 

0.90 

0.22 

0.90 

0.16 

0.8 

0.23 

0.95 

0.38 

0.95 

0.35 

0.9 

0.36 

1.00 

0.47 

1.00 

0.54 

1 .0 

0.54 

1 .05 

0.58 

1 .05 

0.71 

1.1 

0.69 

1.10 

0.76 

1 .10 

0.81 

1.2 

0.76 

1.15 

0.81 

1.15 

0.93 

1 .3 

0.83 

1.20 

0.92 

1 .20 

0.96 

1.4 

0.92 

1 .25 

0.93 

1 .25 

0.99 

1.5 

0.94 

-- 

-- 

1 .6 

0.97 

-- 

-- 

-- 

-- 

1 .7 

0.99 

-- 

-- 

-- 

-- 
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?-54-75-< 


Figure  7.  GRAPHICAL  PRESENTATION  OF  TABLE  3 
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Table  4.  PROBABILITY  OF  WIN  AS  A FUNCTION  OF 
FORCE  RATIO  (Unequal  Breakpoints) 

[q b = 0.67;  qr  = 0.79] 


f 

P ( Red  Wins) 

f 

P(Red  Wins) 

Bo  ' 

R0  = 50 

Bo  ‘ 

Rq  = 100 

1 .00 

0.16 

1 .00 

0.08 

1 .25 

0.24 

1 .25 

0.20 

1 .50 

0.47 

1.50 

0.45 

1.75 

0.60 

1 .75 

0.68 

2.00 

0.73 

2.00 

0.73 

-- 

-- 

2.25 

0.88 

- - 

2.50 

0.94 

Bq  = 50,  Rq  = 75 

Bo  = 

R0  = 200 

0.75 

0.00 

1 .00 

0.00 

0.90 

0.04 

1 .20 

0.10 

1 .05 

0.11 

1 .40 

0.29 

1 .20 

0.22 

1 .45 

0.29 

1 .35 

0.24 

1 .50 

0.40 

1 .50 

0.42 

1 .55 

0.45 

1 .65 

0.58 

1 .60 

0.56 

1 .80 

0.67 

1 .65 

0.58 

1 .95 

0.76 

1 .70 

0.70 

2.10 

0.71 

1.75 

0.71 

-- 

-- 

1 .80 

0.80 

-- 

-- 

2.00 

0.86 

-- 

-- 

2.20 

0.94 

-- 

-- 

2.40 

0.96 

-- 

2.60 

0.99 

-- 

-- 

2.80 

1 .00 
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P (RED  WINS) 
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LANCHESTER  LINEAR  FORCE  RATIO,  f 

2-26-75-5 


Figure  8.  GRAPHICAL  PRESENTATION  OF  TABLE  4 
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qb  = 0.67 


qr  = 0.79  • 

These  figures  are  taken  from  the  ATLAS  model,  where  Red  is  the 
attacking  side  and  Blue  is  the  defending  side. 

The  results  are  qualitatively  the  same  as  for  the  case 
of  equal  breakpoints;  the  reader  may  form  and  interpret  his 
own  conclusions.  Of  interest  is  the  problem  of  solving  the 
equation 

Pl(f)  = 1/2  , i > 0 . 


From  the  data,  it  appears  that  the  solution  is  possibly  inde- 
pendent of  i and  equal  approximately  to  1.57.  Noting  that 


1 - 0.67 
1 - 0.79 


1.57143  , 


we  are  led  to  formulate  the  following  assertion. 


CONJECTURE.  For  all  values  q^  and  q^  of  the  breakpoints 
and  all  values  of  the  force  effectiveness  parameters. 


P 


1/2 


for  all  i. 

Hopefully,  a proof  will  be  devised  soon. 

In  Figures  9,  10,  11,  and  12,  we  give  various  simulated 
probability  distributions  for  numbers  of  survivors  in  these 
situations.  Figure  9 gives  the  distribution  of  the  number  of 
survivors  on  the  winning  side,  for  Rq  = Bq  = 100,  f = 1,  and 
q = q =0.5.  We  note  that  the  narrowest  victories  are  the 
most  likely.  The  expected  surviving  strength  is  approximately 
55  percent.  Figure  10  presents  analogous  data,  but  with 
f = 1.25.  In  this  case,  the  defender  (the  weaker  side)  is 
likely  to  win  only  by  a narrow  margin,  while  the  attacker's 


24 


R = B = 100 
o o 

f = 1 

qb=qr=0-5 

50  55  60  65  70  75  80 


PERCENT  SURVIVOR  STRENGTH  OF  WINNER 


Figure  9.  PROBABILITY  DISTRIBUTION  OF  NUMBERS  OF 
SURVIVORS  (Equal  Forces) 


50  55  60  65  70  75  80 


PERCENT  SURVIVING  STRENGTH  WHEN  VICTOR 


Figure  10.  PROBABILITY  DISTRIBUTION  OF  NUMBERS  OF 
SURVIVORS  (Force  Ratio  = 1.25) 
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Figure  11.  PROBABILITY  DISTRIBUTION  OF  NUMBERS 
OF  SURVIVORS  (Unequal  Breakpoints 
and  Force  Ratio  * 1 ) 


Figure  12.  PROBABILITY  DISTRIBUTION  OF  NUMBERS  OF 
SURVIVORS  (Unequal  Breakpoints  and 
Force  Ratio  = 1.57) 
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most  likely  surviving  strength  is  on  the  order  of  55  percent, 
with  expected  surviving  strength  of  approximately  60  percent. 

The  data  in  Figures  11  and  12  are  based  on  the  ATLAS 
breakpoints  and  on  force  ratios  of  1 and  1.57  (approximate 
force  equality,  in  the  sense  that  each  side  is  equally  likely 
•to  win),  respectively.  The  reader  can  form  his  own  interpre- 
tations and  draw  his  own  conclusions.  We  remark  again  that 
overwhelming  victories  are  not  likely.  In  comparison  with 
data  concerning  the  stochastic  square-law  attrition  process 
presented  in  Reference  [5],  the  margins  of  victory  are  uni- 
formly less  in  the  linear-law  case.  These  margins  are  ex- 
pected, in  the  sense  that  our  intuition  concerning  the  two 
processes  is  that,  in  the  square-law  process,  the  numerically 
superior  side  can  bring  its  full  superiority  to  bear  against 
the  opposition  but  that  it  cannot  do  so  in  the  linear-law 
process . 
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